General plane wave mode functions for scalar-driven cosmology 
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We give a solution for plane wave scalar, vector and tensor mode functions in the presence of any 
homogeneous, isotropic and spatially flat cosmology which is driven by a single, minimally coupled 
scalar. The solution is obtained by rescaling the various mode functions so that they reduce, with a 
suitable scale factor and a suitable time variable, to those of a massless, minimally coupled scalar. 
We then express the general solution in terms of co-moving time and the original scale factor. 
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1. Introduction: On the largest scales the universe 
is homogeneous and isotropic 0. When results from the 
WMAP satellite are combined with other data the energy 
density of the universe as a fraction of the critical density 
is found to be ^Itot — 1-02 ± .02 2]. This is consistent 
with the high degree of spatial flatness predicted by most 
models of inflation 3]. It is therefore safe to assume 
that the invariant element relevant to cosmology takes 
the form, 



ds^ = -dt^ + a\t)dx ■ dx . 



(1) 



The scale factor a{t) is not observable but the ratio of 
its current value to its value at past time t gives the cos- 
mological redshift (plus one) experienced by light emitted 
at that time and received now. 



a{t) 



1 



(2) 



It is common to use redshift as a time variable in cosmol- 
ogy, even for epochs from which we detect no radiation. 
For example, the cosmic microwave radiation was emitted 
within about 100 redshifts before and after Zdec = 1089, 
and the transition from radiation domination to matter 
domination occurred at Zeq ~ 3200 Q • Primordial infla- 
tion is conjectured to have ended at z > e^". 

The scale factor's logarithmic derivative defines the 
Hubble parameter, 



Hit) 



(3) 



Its physical meaning is the rate at which distant matter is 
receding due to the expansion of the universe. Combining 
WMAP with other data sets determines its current value 
to be, Ho = (71 tt) rWT I2- ^he weak energy condition 
implies that the Hubbfe parameter can never increase. 
One typically assumes it has fallen to Hq from an initial 
period of inflation during which it was up to 55 orders of 
magnitude larger. 



The scale factor's second time derivative defines the 
deceleration parameter. 



, ^ ad H 



H^ 

It is less well measured than H{t), but Hubble plots of 
Type la supernovae at high redshift are consistent with 
a current value of go ^ —-6 4]. This phase of nega- 
tive deceleration is a relatively late phenomenon, having 
set in about z ^ 1. Before this it was near the value 
g = +i characteristic of a matter dominated universe. 
At redshifts much larger than z^q ^ 3200 the decelera- 
tion parameter was near the value q — +1 of a perfectly 
radiation dominated universe. And during the conjec- 
tured epoch of primordial inflation [z > e^'^) it would 
have been near the value q = —1. 

The fact that the deceleration parameter has changed 
so much over the course of cosmic history is frustrating 
because many of the phenomena which we can observe 
from the epoch of matter domination {q — +i) are be- 
lieved to have had their origins in quantum fluctuations 
during the epoch of primordial inflation (q ~ —1) j^. 
The various effects can be understood using free quan- 
tum field theory but the mode functions are known only 
for the case of constant q{t). Hence one must employ 
either perturbative or numerical methods to connect the 
inflationary mode functions — whose normalization is 
fixed — with their matter dominated descendants whose 
effects can be observed 0, 111 ■ 

The purpose of this paper is to present an exact solu- 
tion for the scalar, vector and tensor mode functions of 
gravity plus a single, minimally coupled scalar with any 
sort of potential. That we can do this derives from previ- 
ous work in which an exact solution was obtained, for 
any a(t), for the mode functions of a massless, minimally 
coupled scalar, 

+ H\t) [x% k) + ^qit) I] nit, k) = 0. (5) 

The dimensionless parameter x{t, k) = k/aH is the phys- 
ical wave number expressed in Hubble units. We seek a 
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solution (O which is conventionally normalized, 

u{t, k)u*{t, k) — u{t, k)u*{t, k) — i , 



(6) 



and initially (at t ^ ti) positive frequency. The desired 
solution takes the form of a row vector containing the two 
Bessel functions of order = q^^(t), which comprise 
(most of) the solution for constant q{t), multiplied into a 
2x2 matrix that mixes the two solutions whenever q{t) 
evolves, 

u(t,fc) = ' (^-icsc(i/7r)J_^(-|) , J^(-f)) 



xM{t, ti, k) 



l + icot[i^(ti)7r] 



(7) 



The transfer matrix M.{t, ti, k) is the time-ordered prod- 
uct of the exponential of a line integral, 



M{t,U,k) EE exp 



dt'A{t', k) 



(8) 



E 

n=0 



dti I dt2 ■ 

ti Jti 



dtnA{ti,k)---A{tn,k). (9) 



The exponent matrix A{t, k) vanishes whenever q(t) is 
constant. It has the form, 



A{tk)^-v( csc(j^7r)c^(-|) -2id^{-^) 

4 \-2icsc\vT:)K{-^) -csc(i^7r)c,,(-|) 



where the various coefficient functions are 

^2n-2v 



(10) 



1 ^ (-l)"r(n-z.-l)z^"-^-(n-^)-i 
2v^^^ T{n)T{n-v+l)T{n-2v + l) ' 



4 r 
c^{z) — sin(t^7r) i^{v) — 1 — In(iz) 

TT L 



(-l)"r(n- 



-i)z2"n-i 



d^{z) = 



^{n+'y)T{n+l)r{n~iy+iy 
1 (-l)"r(n+;/-^)z'"+^"(n + ;y)-i 



(12) 



,(13) 



Here ipiz) = r'(.z)/r(z). 

Equation ^ is of phenomenological interest because 
it gives the mode functions of linearized gravitons [lO| . 
However, the dominant perturbations imprinted in the 
cosmic microwave background seem to be from the scalar 
modes of gravity plus a minimally coupled scalar 0. The 
perturbative background and the equations defining the 
mode functions of this system are reviewed in Section 2. 
In Section 3 we solve for the mode functions by altering 
the scale factor and time parameter in 0. In Section 
4 the general solutions are re-expressed in terms of the 
actual scale factor and the co- moving time t. We give 
simple and accurate approximations for the ultraviolet 
regime {x{t,k) ^ 1) and the infrared regime {x{t,k) <^ 
1) in Section 5. Section 6 discusses applications. 



2. Gravity with a minimally coupled scalar: The La- 
grangian in which we are interested is, 



C = 



1 



16nG 



^-^d.^d^^g'^'^V^-Viip)^ . (14) 



This system can be solved perturbatively in a variety of 
different field variables a nd gaug e conditions . We will 
use a recent formulation [ill Il2| based upon generalizing 
the de Bonder gauge condition typically employed for 
quantum gravity in flat space. 

The full scalar and metric fields are expressed in terms 
of their background values plus quantum fields. 



f{r],x) = ipo{ri) + 0(77, f) , 



(15) 
(16) 



Here = IGttG is the loop counting parameter of quan- 
tum gravity. As usual, indices on the graviton field ip^^ 
are raised and lowered using the Minkowski metric rj^^. 

Note that we consider the fields to be functions of con- 
formal time rj, 



r] = Vi + 



dt' 



d_ 

drj 



ait) 



dt 



(17) 



This in no way precludes expressing the relations between 
the background fields ¥'o('7) s,iid a{t) in terms of co-coving 
time derivatives, 



-2H-'iH'^ 



^g{^^1 + V{^o)) , (18) 

(19) 



87rG(^-^^-y(^o) 



As usual, an overdot denotes differentiation with respect 
to t whereas a prime stands for differentiation with re- 
spect to 77. 

It is traditional with this system to consider the scalar 
potential V{(f) to be a known function from which the 
scale factor is inferred. However, for our purposes it is 
more convenient to regard a{t) as the known function. 
If one desires the scalar and its potential they can be 
reconstructed using the relations. 



and V{(Pq) 



{2-q)H' 



The generalized de Bonder gauge condition is [lll |. 



(20) 



Ff, = a 



r^,,,~^r.,^-2aH^^o + 2dl^l+qHa^=0 . 

(21) 

In this gauge it can be shown that the linearized fields of 
this system — and the Green's functions needed to solve 
at higher orders — are described by three plane wave 
mode functions Qj{r],k) [l^ . 



Qi^O for / = A,B,C , (22) 



3 



The three parameter functions Oj are, 



1 



3b = — and 
a 



1 



with respect to the co-moving time t. This suggests that 
we define associated Hubble and deceleration parameters, 



7C = 



(23) 



We normalize the various mode functions to obey, 

QiQ}' - QiQ*i = I ■ (24) 

The linearized fields in de Bonder gauge (|21|l can be 
expressed using two mode sums !l2j, 



H =1^ 
9i dti 



and 



-1 



1 dHi 



Hj dti 



A little calculus establishes the identity. 

2 ^r/ d ^2 



\dri 



I 1 +Hi 



\dti/ 



dti 



(33) 



(34) 



1 f d^k 



d^k 
(2^ 



'Y{k) 



(25) 



Applying this to the function 6i itself we can re-express 
the quantity ~6'j /9i which appears in the Qi mode equa- 
tion 



- J^Y.{<^^AkA)QAe''^■'a{k,X)+c.c]. (26) 

Note that in these mode sums we have suppressed the 
arguments of previously defined functions such as a(t), 
H{t), q(t) and Qiiji, k). The scalar and graviton creation 
and annihilation operators are canonically normalized, 



ei 



(35) 



In keeping with the geometrical analogy, we define the 
"physical" wave number in "Hubble" units as. 



XI 



eiHi 



(36) 



Y{k), YHk') = {2T:f5^{k-k') , 



(27) 



The final step is to rescale Qi by a factor of \/6i 
d^ui r 9 3 



a(fc,A), at(A:',A') = {2nf 5^ {k ~ k! ) 5 . (28) 



As usual, graviton polarization tensors are transverse 
{kieij{k,X) — 0), traceless {Sijeij{k, X) — 0), and or- 
thonormal (e^ (fc, A)e*j (fc, A') = S\x')- The other lin- 
earized fields are, 



Ul 



'6iQi 



dt] 



-.11 



ui = 



(37) 

This is the same equation for a plane wave mode function 
ui evolving in the geometry dsj = —dt'j + 9'jdx • as © 
is for the mode function u(t, k) evolving in the geometry 
ds^ = —dt^ + a^dx ■ dx. Even the Wronskian agrees up 
to a sign. 



d 



^Po^{r|.x) = 0, (30) 
ipi-i{r],x) = ^^JtpoQ{■l],x) + tpfJr],x) . (31) 



(^aipao{r],x)j , (29) du} dui 

,„„N dtj dtj 



* I 



QiQ*i — Q'iQ*i 



= ±i 



(38) 



We can absorb the sign by making the B and C modes 
scale to Q*, which does not change ||27J), 



Although the B modes are not required at linearized or- 
der it is no trouble to get them, and they would be ap- 
pear in the retarded Green's functions needed to obtain 
the next order solutions [l^ . 

3. Reduction to the MMC scalar mode functions: Our 
procedure is to define a new scale factor and a new time 
coordinate so that the Qi mode equation (|22|l reduces to 
the mode equation (jSJ of the massless, minimally coupled 
scalar. There does not seem to be any way to motivate 
why this works, but the trick is simple enough that we 
shall simply state it. 

We define the new time variable t/ as follows. 



UAitAjk) = 
UB{tB,k) = 



OcQhiv,k) 



(39) 
(40) 
(41) 



dti = ±6idri 



- = ±9i — 
drj dti 



(32) 



The plus sign is used for functions 9i which grow with 
conventional time (hence I — A) while the minus is taken 
when 9 1 falls (hence / = B,C). Note that 9i plays the 
same role with respect to ±t/ as the scale factor a{t) does 



It follows that ui must be the same function of ti , with 
scale factor 9i, as u{t, k) is of t, with the scale factor a{t). 

4-. Expressing the solution in conventional form: Of 
course physics is not based on the fictitious geometry 
dsj = —dt] + 9jdx ■ dx, so we must express the solutions 
ui as functions of co-moving time t and as functionals of 
the true scale factor a{t). From their definitions (|23() we 
already know how the functions 9i(t) depend upon the 
co-moving time. To determine Hi(t) and qi(t) we need 
the differential relation between ti and t. By comparing 
the transformation H32|l from rj to ti with the transfor- 
mation H17() from 77 to f we infer. 



d 
drj 



dt 



d 
'dTi 



±- 



'9i dt 



(42) 
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As before, the + sign pertains for A modes and the — 
sign for B and C modes. 

Now apply to the parameters of the dsj ge- 
ometry. The A and B mode parameters are simple, 



Ha = H 
Hr = Ha^ 



qA = q , 

IB = na , qb = -q ■ 
The C mode parameters are more complicated 

r 

q , 



He = {l+r)^/l+^Ha^ , qc = 



H 



-r (l + r)2 
where we define the additional parameter, 
1 d 



(43) 
(44) 



(45) 



(46) 



However, note that whenever q{t) is constant — which 
is the case for all familiar cosmologies — the C mode 
parameters degenerate. 



q = He = \/l+qHa^ and qc = —q ■ 

The rescaled wave number is in all cases simple, 

x{t, k) 



xa = a^(i, ^) , ^B — 2;(i, fc) and X(j 



l+r(t) 



(47) 



(48) 



In fact, whenever q(t) is constant, xc = x{t, k). 

We have reached the point where the parameters of 
the dsf geometry can all be considered to be functions 
of co-moving time t: 9i{t), Hi{t), qi{t) and x/(t, k). We 
make the additional definition vjil:,) = ^—q~J^{t). Based 
on the work of the previous section the I ^ A, B,C mode 
functions can be read off from (O by replacing the pa- 
rameters of the ds^ geometry with the corresponding dsj 
parameters, considered as functions of t, 



ui{t, k) \ 



icsc(i^/7r) J_^X-fi) , Juji-^) 
1 



2k V 2qiJ 

xMiit,U,k) , . , , 

As with u(t, fc), the transfer matrix takes the form, 

dt'Ai{t',k) 



Mi{t, ti, k) = P < cxp 



Its exponent matrix obeys the same scheme. 



(49) 



(50) 



TT . / csc(^'/7r)c^,^ff) -2iduj(-^) 



where the coefficient functions (jllllBII are unchanged 
Note that we were able to express the transfer matrix in 
terms of integrations over the co-moving time owing to 
the relation. 



/ dti-^F{xi,qi)= dtviF{ 



XI, qi) 



(52) 



5. The ultraviolet and infrared regimes: The ultraviolet 
regime is defined by x(t,k) ^1. In this limit we know 
from previous work [9| that the three mode functions 
approach the form. 



ui{t, k) 



x-»i V 2k 
where the phase x/ is, 

Xi{t,li) = - 



exp 



iXiiU, k)~ik 



ti{U) 



9i 



Xj{t, k) i — qi{t) TT 



qiit) 



qi{t) 2 



(53) 



(54) 



Recalling dtj — ±9idt/a, and the definition (|17|l of con- 
formal time, we see that the various Q/'s approach, 



Qi{v,k) 



1 



2k 



(55) 



Note that the ± signs in the relation between dtj and 
dt exactly cancels the variations H39I41|I in the relations 
between ui{t, k) and Qi{ri, k). 

The infrared regime {x 1) is more subtle and more 
interesting. We assume that the wave number was ini- 
tially ultraviolet, i.e., x{ti,k) ^ 1. During inflation 
x(t, k) falls, typically exponentially. This drives the 
lower ultraviolet modes through first horizon crossing 
{x{ti, k) = 1) into the infrared regime of x{t, k) <C 1. In 
this regime one of the Bessel functions will enormously 
dominate the other. It is therefore good to re-express the 
general solution in terms of solutions which are the pure 
ibz// Bessel functions at first horizon crossing. 



ui{t, k) = ^ [uj {t, ti,k), uj{t, ti, k)^ 



xA<Kti,^.fc)(i^,,„,[,^(,^),]).(56) 



Here uf{t, ti, k) are the pure J±ui solutions at first hori- 
zon crossing, evolved to time t > ti. 

As long as the condition xi{t,k) ^ 1 is satified, the 
transfer matrix can be worked out up to negligible cor- 
rection factors of order x^ . To the same fractional error 
we can also use the leading terms in the series expan- 
sions of the Bessel functions. The resulting form for the 
— solution is, 



Uj {t, k) 



\/nk 



2 qi' 



(57) 



(51) 

The -|~ solution involves an integral, 



ij{t,k) 



^9f{t) 



2viqiHi9] 



9/ (^r~ 

^ 2 qil 
* dt' 



(58) 
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For the A mode it is the solution that dominates. 
We express UA{t,k) with the standard normaUzation 
times two /c-dependent correction factors, 



UA{t,k) 



ai{t)C,A{k)C,A{k) . (59) 



Hi 



1 



Here Hi = H{ti) is the Hubble parameter at first hori- 
zon crossing. Note that it can depend upon k because 
the time of first horizon crossing depends upon the wave 
number. Indeed, some modes never experience horizon 
crossing. 

The correction factor in that depends upon the 
system's state at first horizon crossing is. 



It is the solution that dominates the infrared C 
mode, 

uc{t, k) 

X I dt'a{t')[l + q{t')]Cic{k)C,c{k). (67) 
The correction factor at horizon crossing takes the form, 

Cic(fc) 



\ e 'ic cos(^)(l + ri) 'ic 



CiAik) 



~ f ^-vii 



(60) 



(68) 

where q\c = qc{ti) and ri = r(ti). We remind the 
reader that the C deceleration parameter and r{t) are, 



Here qi = q{ti) is the deceleration parameter at first 
horizon crossing. The correction factor depending upon 
evolution from to ti is, 

C,A(k) = M]l{ti,U,k)+M'^{ti,t,,k)e'^sec{f-) . (61) 



qc{t) = - 
r{t) 



(69) 



l+r{t) [l + r(t)]2 ' 

^|i„(vTT^)^ (70, 

The other correction factor depends upon qic = qciU), 



Here qi = q{ti) is the initial value of the deceleration Cic{k) = Mf^{ti,ti,k) + M'^{ti,ti,k) e ''-csec(^) 



parameter. Since i'{t) — is typically small during 
inflation, it ought to be a very good approximation to 
simply take the first several terms of the series expansion 
of the transfer matrix in estimating CiAik). Combining 
()59|l and H39|l we see that the infrared limit of Qa takes 
the form, 



(71) 

Taking (|67|l together with H41() gives the last of our in- 
frared limits. 



Qc{v,k) 



Hi 



1 



1 



x<i y/2k{l + qi) y/l+q{t) ait) 



QA{v,k) 



-iH 



^ait)CiAik)C,A{k) . (62) 



For the B mode it is the solution that dominates, 
and only the integral is significant. 



usit, k) 



Hi 1 



dt'a{t')CiB{k)C,B{k). (63) 



\/2kas{t)Jti 
The correction factor from horizon crossing is, 

CiB{k) ^ e f-^ cos(^) . (64) 

The other correction factor depends upon evolution up 
to this time, 

C^B{k) = M^iiti,U,k) + M^iiiti,U,k) e^'^seci f) . 

"(65) 

And we can combine (|63() and (|40|l to extract the infrared 
limit of Qb, 



QsiVik) 



Hi 1 



2k a{t) 



dt'a{t')Cls{k)CUk). (66) 



X dt'a{t')[\ + q{t')]Clc{k)Clc{k). (72) 

6. Applications: The formalism we have developed 
might seem intimidating, and one must of course approx- 
imate the transfer matrix in order to compute anything. 
However, it is often of great value to possess general ex- 
pressions. For example, one can search for unexpected 
effects such as the consequences of the universe having 
passed through q = at the end of inflation. We can 
also explore cosmologies for which the C mode deceler- 
ation parameter H69|) differs appreciably from its small 
r{t) limit of -q{t). 

Our formalism is also useful. For example, our infrared 
limiting form (|72|l allows one to compute the Sachs- Wolfe 
contribution to the primordial scalar power spectrum for 
any evolution a{t) which does not compromise the in- 
frared limit. It should be noted, for non-experts, that this 
primordial power spectrum is almost perfectly flat. The 
rich structure imag ed by WMAP derives from a vari- 
ety of late-time effects which occur after second horizon 
crossing, when the stress-energy is no longer dominated 
by the inflaton and our scalar mode functions are not 
relevant. (However, note that our tensor mode functions 
continue to apply. They depend only upon a(t), not what 
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caused it.) The primordial power spectrum sets one of 
the initial conditions for this late-time cosmology, whose 
dynamics are well understood. 

Assuming the cosmic microwave radiation radiation 
was emitted during pure matter domination {q — +^), 
our result for the primordial scalar power spectrum is 

El, 



Vswik) = ^^\\C^c{k)\mc{k)f 
An 1+qi 



(73) 



The standard result is -^GHi / {1 + qi) so the correction 
factors H68(l and H71|l represent improvements. Because 
different conventions exist in the literature we correspond 
'Psw(^) below to the symbol 5{k) used by Mukhanov, 
Feldman and Brandenberger 0|, to the variable VTi{k) 
used by Liddlc and Lyth and to the quantity ^|(fc) 
used by Lidsey et al. 



Vsw{k) 



25, 



S{k)\\ 
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Vnik) 



225 



AUk) 



(74) 



We also specify how it enters the correlation function be- 
tween temperature fluctuations observed from directions 
ei and 62, 



Tr 



Tr 



n) 

I sw 



2iXofe-(ei— 62) 



(75) 



where Xq = x(to, k) is the physical wave number in cur- 
rent Hubble units. 

A very interesting problem we now have the analytical 
power to tackle is the apparent paradox associated with 
the fact that the scalar power spectrum H73|) diverges 
whenever first horizon crossing occurs during a tempo- 
rary period of de Sitter inflation (i.e., qi = —1). It is 
difficult to understand why the result in this case is not 
of gravitational strength, and one would think it should 
not diverge. This has led to controversy in the past 01 
which it is now possible to resolve. 

Finally, these mode functions can be used to study 
back-reaction in the gravity plus scalar system. It has 



been established that no secular back-reaction occurs at 
one loop for simple models 0, . This formalism al- 
lows one to study an arbitrary single-scalar model of in- 
flation. 
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